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Abstract 



We investigate the breakdown of normal hyperbolicity of a manifold of equilibria 
of a flow. In contrast to classical bifurcation theory we assume the absence of any 
flow-invariant foliation at the singularity transverse to the manifold of equilibria. 
We call this setting bifurcation without parameters. 

In the present paper we provide a description of general systems with a manifold of 
equilibria of codimension one as a first step towards a classification of bifurcations 
without parameters. This is done by relating the problem to singularity theory of 
maps. 
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1 Introduction 



We study dynamical systems with manifolds of equilibria near points at which normal 
hyperbolicity of these manifolds is violated. Manifolds of equilibria appear frequently in 
classical bifurcation theory by continuation of a trivial equilibrium. Here, however, we 
are interested in manifolds of equilibria which are not caused by additional parameters. 
In fact we require the absence of any flow-invariant foliation transverse to the manifold 
of equilibria at the singularity. We therefore call the emerging theory bifurcation without 
parameters. 

Albeit the apparent degeneracy of our setting (of infinite codimension in the space of 
all smooth vectorfields) there is a surprisingly rich and diverse set of applications ranging 
from networks of coupled oscillators |Lie97] , viscous and inviscid profiles of stiff hyperbolic 
balance laws |HL05j . standing waves in fluids |AFL08t AFL09J, binary oscillations in 
numerical discretizations |FLA00b] . population dynamics |Far84] . cosmological models 
|WE05j . and many more. The present paper is a first step towards a classification of 
bifurcations without parameters. 

Consider a vector field 

x = f(x,y) GR", 
y = g(x,y) € R m 

with a manifold of equilibria { (0, y) : y G R m }, i.e. 

f(0,y) = 0, g(0,y) = 0. (1.2) 

As long as the manifold remains normally hyperbolic, i.e. the linearization of / on the 
manifold has no purely imaginary eigenvalues, 

specd x f{0,y) n iR = 0, (1.3) 

there exists a local flow-invariant foliation with leaves homeomorphic to a standard saddle, 
for example by the theorem of Shoshitaishvili [Sho75j. Bifurcations are characterized by 
a non-hyperbolic block A of the linearization 



A(y) \ / d x f d y f 
B(y) \d x g d y g 



(0,1/), (1.4) 
say at the origin, i.e. the spectrum of A(0) intersects the imaginary axis, 



specA(O) n iR ^ 0. 



(1.5) 
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Restricting to a center manifold we can assume 

specA(O) C iR. (1.6) 
We will always assume that the vector field is smooth enough to allow suitable expansions. 
Note the analogy to classical bifurcation theory where y would be a parameter, i.e. 



g = 0. For references on classical bifurcation theory see for example |Arn83llHK91[[Kuz95[ 
IVan89j and the references there. In the classical case, g = 0, however, the flow invariant 
transverse foliation {y = const.} is also present in a neighborhood of the bifurcation point. 



This is no longer true in the general case (1.1, 1.2) without parameters. Indeed, generic 



functions g of the form (1.2) yield a drift in the "parameter" direction y which excludes 



any flow-invariant foliation transverse to the manifold of equilibria near a singularity (1.5). 
Thus, the resulting nonlinear local dynamics differ considerably from classical bifurcation 
scenarios. 



A rigorous analysis of bifurcations without parameters (1.1, 1.2, 1.6) has been carried 
out for the following cases: 



1, with linearization 



(i) A simple eigenvalue zero of A (transcritical point), n = m 
at the bifurcation point: 

/ 

V 1 
see |Lie97j . 

(ii) A pair of purely imaginary nonzero eigenvalues of A (Andronov-Hopf point), n = 2, 
m = 1, with linearization at the bifurcation point: 

/ 1 \ 

-10 0, 

V o oo/ 

see |FLA00a] . A partial description can also be found in jFar84j . 

(iii) An algebraically double and geometrically simple eigenvalue zero of A (Bogdanov- 
Takens point), n = m = 2, with linearization at the bifurcation point: 



( 








\ 


1 














1 
















o J 



see |FL01j . With additional symmetries, this case is also studied in |AFL08l [AFL 09J 



Codim-one manifolds of equilibria and bifurcation without parameters 



3 



Note the nonzero blocks B at transcritical and Bogdanov-Takens points yielding a drift 
in y-direction and excluding any flow-invariant transverse foliation. At Andronov-Hopf 
points this drift is induced by a generic second-order term U y A x (^)(0) ^ which is the 
leading order term of the drift in ^/-direction averaged over the periodic motion of the 
linearization. 

Bifurcations without parameters can also appear combined with additional parame- 
ters, for example g(y) = #(2/1,2/2) = {gx{yii 2/2), 0). For Bogdanov-Takens points the case 
of a generic vector field with two-dimensional equilibrium manifold turns out to be equiv- 
alent to the case of a generic one-parameter family of vector fields with one-dimensional 
equilibrium manifolds, at least up to leading order of the suitably rescaled normal form; 
see |FL01j . Both viewpoints are closely related in this case. In the example of a trans- 
critical point with drift singularity, studied in section |2j however, both settings lead to 
drastically different dynamical systems, see Remark |2.2| 

In the present paper we analyze the case x G R, y G R m of dynamical systems with 
a codimension-one manifold of equilibria. As it turns out, the dynamics near bifurcation 
points of codimension m on these manifolds can be related to singularities of smooth 
maps h : R m+1 — y R; the Theorem 3^ This correspondence permits the application of 
singularity theory and most notably of the classification of singularities to bifurcations 
without parameters. It might serve as a first step towards a classification of general 
bifurcations without parameters. 

In section [2] we will start with an example to illustrate the general theorem formulated 
and proved in section |3j In section [4] we conclude with a discussion of further questions, 
possibilities and problems. 

Acknowledgement: This work was supported by the Deutsche Forschungsgemein- 
schaft. 



2 Transcritical points with drift singularity 

Let us first discuss the cases m = 1,2 of bifurcations along one- and two-dimensional 
equilibrium manifolds as an example to illustrate the general theorem in section [3j 



4 



Stefan Liebscher 



x 



A 

trivial equilibria 




Figure 2.1: Transcritical bifurcation point: classical (a) and without parameters (b). 
2.1 Transcritical point 

The case m — 1 of a line of equilibria in a 2-dimensional center manifold has already 
been studied in |Lie97] . see also |FL02j . In classical bifurcation theory, the only robust 
bifurcation of 

x = f(x,X) 6R, /(0,A) = 0, 
A = eft 

is the transcritical bifurcation with normal form 

x 



(2.1) 



x{x - A) + h.o.t. (2.2) 

see figure 2.1[ a). It is caused by the nontrivial eigenvalue of the linearization at the 
equilibrium x = crossing zero with nonvanishing speed as the parameter A increases. 
Together with the non-degeneracy condition d^f(0, 0) ^ 0, this implies the above normal 



form (2.2). 



(2.3) 



Without parameters, 

x = f(x,y) GE, /(0 j2 /) = 0, 

V = g(x,y) G E, g(0,y) = 0, 

the nontrivial eigenvalue d x f(0, y) can change sign along lines of equilibria {y = 0}. The 
generic normal form reads 

x = xy + h.o.t., 

y = x, 



(2.4) 
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see figure pTTj fo). It requires the same transversality condition of the nontrivial eigenvalue 
as the classical transcritical bifurcation. The non-degeneracy condition, however, is re- 
placed by d x g(0, 0) 7^ and yields the two-dimensional Jordan block of the linearization 
at the transcritical point. Trajectories form parabolas with tangency to the line of equi- 
libria at the transcritical point. The flow direction is reversed on opposite sides of the 
equilibrium line. 

2.2 Parameter dependent transcritical point with drift singu- 



Along two-dimensional equilibrium manifolds we expect transcritical points to form one- 
dimensional curves, by implicit function theorem. At isolated points one of the non- 
degeneracy conditions may fail and codimension-two singularities appear. We will discuss 
the case of failing drift condition, first in a one-parameter-family on lines of equilibria and 
then along a two-dimensional equilibrium surface. 

With one parameter, the setting is as follows. We consider a system 



with the following properties: 

(i) For all parameter values, there exists a line of equilibria, F(0, y, A) = 0, forming a 
plane of equilibria in the extended phase space. 

(ii) For all parameter values, the origin is a transcritical point, i.e. the origin has an 
eigenvalue zero in transverse direction to the equilibrium plane, d x f(0, 0, A) = 0. 

(iii) For all parameter values, this nontrivial eigenvalue crosses zero with nonvanishing 
speed as y increases, d y d x f(0, 0, 0) > 0. 

(iv) At A = the drift non-degeneracy condition fails, d x g(0, 0, 0) = 0. 

(v) This drift degeneracy is transverse, i.e. the drift changes direction with nonvanishing 
speed, as A increases, d\d x g(0, 0, 0) > 0. 



larity 





(2.5) 



A 



The first condition is our structural assumption, (iii,v) are non-degeneracy assumptions 
fulfilled generically, and (ii,iv) describe our bifurcation point. This setup is robust, i.e. 
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under small perturbations of F respecting (i) there is a point near the origin satisfying 
(ii-v) for the perturbed system. From the viewpoint of singularity theory, (ii,iv) define 
a singularity of codimension two that is unfolded versally by the coordinate y along the 
line of trivial equilibria and the parameter A. 

Condition (i) allows us to factor out x, 

F{x,y,X) = xF(x,y,X), (2.6) 
with smooth F. Conditions (ii-v) yield an expansion 

F(x iy ,\) = ( aX + by ) + 0((\x\ + \y\ + \X\n (2-7) 
\ cx + ay + a A I 

with coefficients a, b, c, d, a G R, b > 0, a > 0. We assume an additional non-degeneracy 
condition 



(vi) The matrix 



<W(]f) (o,o,o)= r \ 



is hyperbolic, i.e. has no purely imaginary eigenvalues. 
Setting 

5 := ad — be, r := a + d, (2. 
for determinant and trace, we therefore have 5^0, and r ^ if 5 > 0. 



Applying the multiplier x _1 to system (2.6) preserves trajectories for x ^ but 
reverses their direction for x < 0. After the coordinate transformation x = x,y = ax + by, 
A = bor\, we obtain 

.-^ \A +0((M + M + |A|) 2 ). (2.9) 
—ox + ry + A / 

This yields a bifurcating equilibrium at (x,y) ~ (A/5, 0). Transversality of the branch 
of equilibria with respect to the trivial line of equilibria as well as the hyperbolicity of 
the nontrivial equilibria is ensured by condition (vi). Therefore, terms of higher order in 




(2.9) will preserve this structure. See figure 2.2 for phase portraits in various cases. Note 



the appearance of the generic transcritical bifurcation without parameters, figure 2.1 for 
A ^ 0. 
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bifurcating saddle 




bifurcating focus 




bifurcating node 
Stable set of the origin in green, unstable set in red. 
Figure 2.2: Drift singularity along a one-parameter family of transcritical points. 



8 



Stefan Liebscher 



2.3 Transcritical point with drift singularity without parame- 



Replacing the parameter A discussed above by an additional direction of a plane of equi- 
libria, the drift along this manifold of equilibria is now a two-dimensional vector. If will 
not vanish along generic one-dimensional curves. The drift singularity along curves of 
transcritical points is therefore not characterized by a vanishing drift but rather by a drift 
direction orthogonal to the curve of transcritical points. (A drift in A-direction was not 
possible before.) 

The correct setup is given by a system 



U — (2/I5 2/2) , 9 = {91,92), with the following properties: 

(i) The y-plane consists of equilibria, F(0,y) = 0. 

(ii) There is a transcritical point at the origin, i.e. the y-plane loses normal hyperbolicity 
at this point, d x f(0, 0, 0) = 0. 

(hi) This loss of normal hyperbolicity is caused by the transverse eigenvalue crossing zero 
transversally, V y d x f(0, 0, 0) 7^ 0. Without loss of generality, the gradient points 
in yi-direction, i.e. d yi d x f (0,0,0) > 0, d y2 d x f(0, 0, 0) = 0. By implicit function 
theorem, this gives rise to a curve of transcritical points tangential to the y 2 -axis. 

(iv) At the origin, the drift non- degeneracy transverse to the curve of transcritical points 
fails, 9 x ^i (0,0,0) =0. 

(v) This drift degeneracy is transverse, i.e. the drift direction crosses the tangent to the 
curve of transcritical points with nonvanishing speed along the curve of transcritical 



(vi) The drift does not vanish at the origin, i.e. there is a component tangential to the 
curve of transcritical points, d x g2(0, 0, 0) > 0. 

Note that conditions (i-v) correspond to the conditions of the previous section. Again, 
the degeneracies (ii,iv) are robust under perturbations satisfying (i), provided the non- 
degeneracy conditions (iii,v,vi) hold. The signs of d yi d x f(0, 0, 0) and 9^(0,0,0) in (iii,v) 
can be inverted by reflecting y\ 1— > —y\ and y 2 h-> —yi, respectively. 



ters 





(2.10) 



points, 9^/(0,0,0) 9^9^(0,0,0) + cL 2 2 d x f(0, 0,0) 9^(0,0,0) ^ 0. 



Codim-one manifolds of equilibria and bifurcation without parameters 



9 



The non-degeneracy condition (vi) indeed yields 
d 



V(in,V2) 9 *f> d x9) (0, ^(ya), Jte) 



(2.11) 



3/2=0 



where (x,yi,y 2 ) = (0, 1^(2/2)5 2/2)5 $(0) = 0, $'(0) = 0, is the curve 7 of transcritical points. 
Locally, we could reparametrize y to achieve 1? = 0. Conditions (ii,iii,vi) would then read: 
d x f(Q, 0, y 2 ) = 0, d yi d x f(Q,Q,y 2 ) > 0, d y2 d x gi(0, 0, 0) 7^ 0. But let us continue with the 
general setup. 



As in the parameter-dependent case ( |2.6[ ), we can factor out x due to condition (i), 

(2.12) 



F(x,y) = xF(x,y) = x 



f(?,y) 



However, this time, due to non-degeneracy (vi) no equilibrium remains, 

F(0,0,0) = (0,0,^ 2 (0,0,0)) ^ 0. (2.13) 
Applying the flow-box theorem, there exists a local smooth diffeomorphism 

h(zo,zi,z 2 ) = $ Z2 (z ,z h O), (2.14) 

where <3>t denotes the flow generated by the vector field F. This diffeomorphism fixes the 
origin and transforms F into the constant vectorfield, 

[Dh(z ,z 1 ,z 2 )]- l F(h{z ,z 1 ,z 2 )) = . (2.15) 

V 1 / 

Applying the same transformation to the original vector field F, we obtain 

( ° ^ 

[D/i(^)]- 1 F(/i(z)) = [Dh(z)]- 1 h (z)F(h(z)) = , (2.16) 

V h (z) ) 

where h = (ho, hi, h 2 ). 

In a suitable neighborhood of the origin, the vectorfield F is flow-equivalent to a 
vectorfield 

z 2 = h (z ,zi,z 2 ) (2.17) 

on the real line depending on two (classical) parameters (zq,Zi). Expansion of ho using 
(2.14) and conditions (ii-vi) yields 



z 2 = az 2 + (cqZq + cizi)zl + (bzi + c 2 z + c 3 Zq + c^z z 1 + c 5 zl)z 2 + z + 0(\z\ 4 ) (2.18) 
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Figure 2.3: Cusp singularity z 2 = az 2 + z\z 2 + Zq with a = — 1. Reverse direction of 
trajectories and signs of zq, z\ for a = + 1. The fold line 7 is connected by heteroclinic 
orbits to the curve a, both curves have a common tangent at the origin. 



with 



a 
b 



(d yi dj{0)d y2 d xgi {0)+dld x f{0)d x g 2 {0))d x g 2 (0) ^ 0, 
d yi d x f(0) 0. 



(2.19) 



In particular ho(0, 0, z 2 ) 



az> 



+ C(|z 2 | 4 )- This is a cusp singularity. See |GG73l lGib79[ 



IArn94[ |AGZV83 IMur03j for a background on singularity theory and its connection to 



dynamical systems. In fact, non-degeneracies (2.19) allow to diffeomorphically transform 



(2.18) into the normal form 



z 2 = ±zl + Zl z 2 + z + O(z?), (2.20) 
for arbitrary normal-form order N, see for example [BG92], proposition 6.10. This is a 



minimal versal unfolding of the cusp singularity. See figure 2.3 



Reverting the flow-box transformation, the cusp singularity yields a description of 
the local dynamics near a transcritical point with drift singularity on a two-dimensional 
manifold of equilibria. Note in particular the cusp-shaped fold line 
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27 

T 



± 2 z + O(zo 



JV/3\ 



of the manifold of equilibria that is connected by heteroclinic orbits to the curve 



a 
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Figure 2.4: Transcritical point with drift singularity on a plane of equilibria. Stable set 
of the line 7 of transcritical points in green, unstable set in red, selected trajectories in 
blue. Two different views for a = — 1. Reverse direction of trajectories and switch colors 
of manifolds for a = +1. 



Proposition 2.1 Under condition (i-vi) the vector field (2.10) in a local neighborhood U 
of the origin is flow- equivalent to the cusp singularity ( 2.20\ ). Depending on the sign of 
the cubic term, that is the sign of a = sign (d yi d x f(0) d y2 d x gi(0) + d y2 d x f(0) d x g 2 (0)^j, all 
trajectories in U converge to an equilibrium {0,y) in forward time (a = — 1) or backward 
time (a = +1 ). 

In U , the transcritical points on the manifold of equilibria form a curve 7 through 
the origin. The unstable (for a — — 1) and stable (for a — +1) sets, respectively, of the 
two components 71,72 0/7 \ {0} form manifolds of heteroclinic orbits on opposite sides 
of the manifold of equilibria. Their targets in forward time (a = —1) or backward time 
(a = +1) again form curves o"i j2 on the manifold of equilibria with a\ U {0} U a 2 being a 



tangential curve to 7. See figure 2.4 



Remark 2.2 In contrast to the parameter- dependent drift singularity no equilibria bifur- 
cate. In fact, the drift non-degeneracy excludes any kind of recurrent or stationary orbits 
except the primary manifold of equilibria. 
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3 General Bifurcation at codimension-one manifolds 

We consider the general case of a manifold of equilibria of codimension one, 

i ) = F(x,y) = I V x G R, y G IT. (3.1) 

Typically, such a system will reduced system on a center manifold of finite 

smoothness. Following the discussion in the previous section we obtain the following 
theorem. 



Theorem 3.1 The exists a generic subset of the class of all smooth vector fields (3.1) 
with an equilibrium manifold {x = 0} of codimension one. For every vector field in that 
class the following holds true: 

At every point (x = 0, y) the vector field is locally flow equivalent to a m-parameter 
family 

z m = ±z e + 1 + Y / z k z k m + 0(z%), (3.2) 

k=0 

< £ < m, of vector fields on the real line. Here N is the arbitrary but finite normal- 
form order bounded by the smoothness of the initial vector field (3.1), f,g£ C M , N < M, 



N < oo. This is a versal unfolding of the singularity z m = iz^ 1 at the origin. 

In particular, near bifurcation points of codimension m, that appear robustly at iso- 
lated points on the equilibrium manifold, the vector field is locally flow equivalent to 



m—l 



Z m — i 2 m + + Z k z m + 0( Z m)-i (3-3) 

k=0 

i.e. a universal unfolding of the singularity z m = ±z™ +1 at the origin. 

Proof. The equilibrium condition f(0,y) = g(0,y) = for all y G R m allows us to factor 
out x. 

F(x,y) = xF(x,y) = x\ f }*' V \ I . (3.4) 

The resulting vector field F : R m+1 — > R m+1 does not vanish on the m-dimensional 
submanifold {x = 0}, for generic F. Without loss of generality, consider a neighborhood 
U C R m+1 of the origin. 

We can apply the flow-box theorem to F: Take a local smooth section 

S : R m D V — >U, (3.5) 
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through the origin, S(0) = 0, transverse to the vectorfield F in U. Let Q t be the flow 
generated by F. Then the flow-box transformation 

h(z , ...,z m ) = $ 4ra (E(zb,-,2m-i)) (3-6) 

transforms F into the constant vector field [Dh]~ 1 (Foh) = (0, . . . , 0, 1). Again, $j denotes 
the flow to the vector field F. Applying the transformation h to the vector field F\u, we 
obtain a m-parameter family [D/i] _1 (F oh) = (0, . . . , 0, 7r x h) of vector fields on the real 
line in a neighborhood V of the origin. 

Classification of germs of vector fields and their versal unfoldings is the topic of 
singularity or catastrophe theory. 

Singularities on the real line have the form z m = rkz^ 1 . In generic m-parameter 
families at most m + 1 leading coefficients of the Tailor expansion vanish, i.e. I < m and 

i-i 

Z m = ± Z m X + X! Ck(z , Z m _i)Z k m + 0(z e i ^ 2 ). 
k=0 

The coefficient Q vanishes by linear transformation of z m and the map (zo, z m -i) t— > 
(( , has full rank, generically. Remainder terms, 0(z^ 2 ), can be pushed to any 



finite normal-form order, by a suitable coordinate change. This yields system (3.2). See 
also jBG92j . chapter 6. 

Genericity conditions amount to algebraic conditions of the coefficients of the Taylor 



expansion at the origin. These conditions translate via (3.6) to generic conditions on F. 



The versal unfolding (3.2), one the other hand, is a system of the form (3.1). There- 
fore, it represents the versal unfolding of a generic singularity along m-dimensional man- 
ifolds of equilibria in (m + l)-dimensional phase space. ixi 



4 Discussion 

The present result is a first step towards a more systematic treatment of bifurcations 
without parameters than done by case studies in |FLA00a| IFLOlt IAFL08] . 

The removal of the manifold of equilibria by a scalar, albeit singular, multiplier 
greatly facilitates the analysis but restricts it to the case of manifolds of codimension one 



in the phase space, see (2.12) and (3.4). Hopf points and Bogdanov-Takens points require 
manifolds of equilibria of at least codimension two. Their analysis in |FLA00allFL01j uses 
a blow-up or rescaling procedure reminiscent of the scalar multiplier used here. It seems 
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worthwile to closer connect these bifurcations without parameters to singularity theory. 
This might provide a suitable setting to also include singularities of the set of equilibria 
and generalize the manifold to varieties. 

Contrary to classical bifurcation theory, no recurrent dynamics has been found so 
far near bifurcation points without parameters. For the codimension-one manifolds of 
equilibria discussed here, the drift nondegeneracy yielding the flow-box transformation 
prevents any recurrent dynamics. Similar drift conditions hold true at generic Hopf and 
Bogdanov-Takens points. In fact this is the drift which distinguishes bifurcations without 
parameters from classical bifurcations by preventing any flow-invariant transverse folia- 
tion. Recurrent dynamics should be possible at bifurcations points of higher codimension 
as the drift condition gets less restrictive. 
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